Soft-mode behavior of electromagnons in multiferroic manganite 
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The behavior of the low-frequency electromagnon in multiferroic DyMnOs has been investigated 
in external magnetic fields and in a magnetically ordered state. Significant softening of the electro- 
magnon frequency is observed for external magnetic fields parallel to the a-axis (i3||a), revealing a 
number of similarities to a classical soft mode behavior known for ferroelectric phase transitions. 
The softening of the electromagnon yields an increase of the static dielectric permittivity which 
follows a similar dependence as predicted by the Lyddane-Sachs- Teller relation. Within the geom- 
etry _B1|6 the increase of the electromagnon intensity does not correspond to the softening of the 
eigenfrequency. In this case the increase of the static dielectric permittivity seem to be governed by 
the motion of the domain walls. 



INTRODUCTION 

Soft modes have long been studied in connection with 
structural phase transitions in crystalline solids [l|-l3|. 
Close to such phase transition the frequency of the nor- 
mally lowest lattice vibration strongly decreases '4] (soft- 
ens) as a result of a "flattening" of an effective potential. 
Therefore, soft modes reflect an instability of the crystal 
lattice [5| and are often accompanied by a substantial 
nonlinearity of the system. Especially in ferroelectric 
crystals the softening of a lattice vibration is followed 
by a divergence of the static dielectric permittivity. In 
a simple case a lattice vibration can be written in the 
Lorentz form, 
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Here w^ is the resonance frequency, ujlo is the frequency 
of the longitudinal phonon, 7 is the damping, Soo is the 
high frequency permittivity, and As is the dielectric con- 
tribution of the mode to the static dielectric permittivity. 
From this equation the well known Lyddane-Sachs- Teller 
(LST) relation immediately follows [3| 



£(0)/eoo = (£00 + Ae)/eao = [i^Lo/^r] 



(2) 



Close to structural phase transition in ferroelectrics uJr 
softens and ujlo basically remains constant [l|, which 
leads to the mentioned divergence of the static dielectric 
permittivity. 

Physically similar conclusion can be drawn taking into 
account that the spectral weight of a single excitation is 
often conserved. In spectroscopy this is generally termed 
as a sum rule and can be derived from first principles [6| . 
The spectral weight is proportional to the electron den- 



sity and can be written as S' ~ Ae • 



Evidently, in 



order to preserve the constant spectral weight with de- 
creasing resonance frequency, the dielectric contribution 



must diverge as Ae ~ l/o;^. The last formula works well 
for the soft lattice vibrations and represents another form 
of the LST relation. 

Although governed by different mechanisms, qualita- 
tively similar effects can be observed in glassy materials 
as well. In this case a broad structural relaxation mode, 
called alpha-relaxation [7.] , moves towards zero frequency 
reflecting the freezing of the ionic movement. Contrary 
to the crystalline solids, in glasses the situation is more 
complicated and the spectral weight of the structural re- 
laxation is not conserved. However, the divergence of the 
static dielectric constant close to glass transition has been 
put into discussion for non-crystalline solids as well [8s|. 

Another sum rule can be derived in case of dielectric 
spectroscopy, which of course correlates with the con- 
servation of the spectral weight. This sum rule states 
that the static dielectric permittivity can be calculated 
as a sum of all contributions from absorption processes 
at finite frequencies |6|. Therefore, in experiment one 
might try to always correlate the changes in static prop- 
erties with changes in the absorption spectra at finite 
frequencies. As a characteristic example, the suppres- 
sion of an electro-active mode in multiferroic GdMnOa 
and TbMnOs has been put forward as a spectroscopic 
explanation of field induced changes in the static dielec- 
tric constant [9,]. 

Multiferroics represent an intriguing class of materi- 
als in which electric and magnetic orders coexist |1CH12| . 
Most interesting effects occur if both orders are strongly 
coupled. This leads to such effects like the control of elec- 
tric polarization by external magnetic field or the control 
of magnetization by electric field. In addition to static 
polarization, dynamic properties of multiferroics are very 
rich and in various compounds show the existence of a 
series of new excitations. These excitations are electri- 
cally active magnetic modes of the cycloidal spin struc- 
ture and they have been called electromagnons 13|, [iJ] . 



The electromagnons may be suppressed in external mag- 
netic field leading to substantial changes of the dielectric 
permittivity in a broad frequency range. 

DyMnOa belongs to the most studied multiferroic 
manganites with orthorhombic structure. Below the Neel 
phase transition at T = 39 K DyMnO^first possesses an 
incommensurate magnetic order |15l4l7| . For T < 19K 
this order turns into a spin cycloid with the manganese 
spins rotating in the crystallographic be plane. As has 
been proven both theoretically and experimentally, this 
spin structure leads to the occurrence of the static electric 
polarization parallel to the c axis 16|, |l8J, [iQj . Similar to 
such multifcrroics like GdMnOs or TbMnOa, DyMnOa 
shows the series of electromagnons at finite frequencies, 
which basically consists of two modes at 2 and 6 meV (15 
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and 50 cm~^), respectively IJ, l20|, |2l|. Although the 
physical mechanism of the electromagnons is not fully 
understood until now, the high frequency mode seems to 
correspond to a zone boundary magnon 22|, [23[. This 
magnon acquires electric dipole activity and becomes 
visible in the optical spectra as a result of the Heisen- 
berg exchange mechanism combined with the cycloidal 
spin structure. The up-to-date situation with the low- 
frequency electromagnon is a bit more complicated. Ac- 
cording to recent experimental results [2J, |25| on the 
closely similar multiferroic TbMnOa, the low- frequency 
electromagnon corresponds to an eigenmode of the cy- 
cloidal spin structure 26] which becomes infrared active 
due to an incommensurate spin modulation of the cy- 
cloid. We note that an alternative explanation based on 
anisotropic effects has been suggested as well [27|, l28| . 

In this work we present the comparison of the electro- 
magnon dynamics in external magnetic fields with the 
changes in the static dielectric permittivity in the mul- 
tiferroic DyMnOa. Substantial decrease of the electro- 
magnon frequency is observed in the cycloidal magnetic 
phase. This softening behavior correlates with an in- 
crease of the static dielectric permittivity, thus revealing 
similar dependence as predicted by the Lyddane- Sachs- 
Teller relation. 
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FIG. 1: Frequency dependence of the real (a,b) and imagi- 
nary (c,d) parts of the dielectric permittivity e* = £i -I- ie2 
of DyMnOa along the a axis for different temperatures and 
external magnetic fields. For increasing magnetic fields a soft- 
ening and a growth of intensity of the electromagnon is clearly 
observed. The inset shows the (B,T) phase diagram for static 
magnetic fields along the a axis. IC - incommensurate mag- 
netic state, "afc" and "6c" denote the a6-plane and the 6c- 
plane oriented cycloids, respectively. Error bars for the bc- 
ab transition reflect the hysteresis and history dependence in 
different experiments, e indicates the ac electric field of the 
electromagnetic wave. 



trol alpha-analyzer) was used for static dielectric mea- 
surements. Single crystals of DyMnOs have been grown 
using the floating-zone method with radiation heating. 
The samples were characterized using X-ray, magnetic, 
dielectric and optical measurements 13|, |3l| . The results 
of these experiments including the magnetic phase dia- 
grams are closely similar to the published results 16 1. 



RESULTS AND DISCUSSION 
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EXPERIMENTAL DETAILS 

The spectroscopic experiments in the terahertz fre- 
quency range (3 cm~^ < t/ < 40 cm~^) have been 
carried out in a Mach-Zehnder interferometer arrange- 
ment 29l |30| which allows measurements of the ampli- 
tude and the phase shift in a geometry with controlled 
polarization of radiation. Dynamic dielectric properties 
e*{T,B) = £i + 162 were calculated from these quan- 
tities using the Fresnel optical equations for the com- 
plex transmission coefficient. The experiments in ex- 
ternal magnetic fields up to 8 T have been performed 
in a superconducting split-coil magnet with polypropy- 
lene windows. A frequency-response analyzer (Novocon- 



Figure [T] shows typical spectra of the low frequency 
electromagnon in DyMnOs, which have been obtained 
in the magnetically ordered state with the spin cycloid 
oriented in the crystallographic be plane. This phase is 
indicated as "6c" in the phase diagram shown in the inset 
to Fig. [1] Application of an external magnetic field par- 
allel to the crystallographic a axis leads to a rotation of 
the be cycloid around the b axis from the be plane to the 
ab plane. The spectra in Fig. [T] clearly demonstrate that 
in external magnetic fields B\\a the electromagnon shifts 
to lower frequencies and gains intensity. This behavior 
reveals already at this point a close similarity to classi- 
cal soft modes. In order to investigate this similarity in 
more details, we have carried out the qualitative analysis 
of the electromagnon using the Lorentz-oscillator model 
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FIG. 2: Comparison of the static and dynamic properties of 
DyMnOa along the a a:xis. Diamonds represent the measured 
dielectric contribution of the electromagnon. Circles - dielec- 
tric contribution of the electromagnon as predicted by the 
LST relation. Triangles - static dielectric permittivity with a 
high-frequency value £oo ~ 25 subtracted. 



and compared these results with the static dielectric per- 
mittivity. 

Figure [5] shows the magnetic field dependence of the 
dielectric contribution of the electromagnon in compari- 
son with the static permittivity. Blue diamonds represent 
the field dependence of the dielectric contribution Ae of 
the low-frequency electromagnon showing an increase by 
more than a factor of two as approaching the phase tran- 
sition to the a6-oriented spiral around 7 T. The dielectric 
contribution closely correlates with the decrease of the 
resonance frequency, which is demonstrated by plotting 
Ae • ^\Xi (red circles). This plot corresponds directly to 
the Lyddane- Sachs- Teller relation and reflects the con- 
servation of the spectral weight of the electromagnon in 
external fields parallel to the a axis. Open triangles in 
Fig. [2] show the field dependence of the static dielectric 
constant in DyMnOa as measured at 10 kHz. Here we 
subtracted the contribution from the higher frequency 
processes (Eoo = 25) which is given by electronic tran- 
sition, phonons, and a second electromagnon [20| at 6 
meV. We observe a close correlation between static and 
dynamic properties in spite of more than seven orders of 
magnitude difference in frequency. According to the sum 
rules mentioned above, this result demonstrates that for 
the geometry i3||a the changes in the static properties 
are nearly completely governed by the softening of the 
electromagnon and no other contributions exist between 
kHz and THz frequencies. 

We note that a correspondence of lattice soft modes 
and electromagnons is not straightforward. The latter 
are defined in the ordered state and reveal critical be- 
havior approaching a phase transition as function of mag- 
netic field. The LST-derived equation for DyMnOs re- 
lates the dielectric constant at zero magnetic fields to 
that in finite fields, while it relates transverse and longi- 
tudinal modes in dielectrics. The softening of the elec- 
tromagnon in external fields can be qualitatively under- 
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FIG. 3: Frequency dependence of the real (a,b) and imaginary 
(c,d) parts of the dielectric permittivity of DyMnOs for e|ja 
for different external magnetic fields i3||fe and at T — 6K. 
The inset shows the (B,T) phase diagram for static magnetic 
fields along the b axis. Notations are the same as in Fig. [l] 



stood taking into account the switching of the orientation 
of the spin cycloid. Similar to many other structural 
transitions the effective stiffness of the cycloid probably 
tends to zero on the phase border between the be- and 
ab cycloids, leading to a softening of the electromagnon. 
The unresolved question is: why the spectral weight of 
the electromagnon is conserved during the softening of 
the eigenfrequency? In case of classical softening of the 
lattice vibration one normally argues that the spectral 
weight of the soft mode is directly connected to the to- 
tal number of electrons in the material. In agreement 
with the charge conservation a constant spectral weight 
may be expected for soft phonons. In case of a magnetic 
cycloid the electromagnon gains the spectral weight as 
a result of a complex interplay of various mechanisms. 
Therefore, we cannot use the conservation of the mag- 
netic moment as an argument, and the observed conser- 
vation of the spectral weight remains an actual problem. 



B\\b 



We turn now to the experimental geometry in which 
the transition from the be- to the ab cycloid is achieved 
by magnetic fields along the b axis. Although here some 
hints to a soft mode behavior could be detected as well, 
the results turned out to be more complicated to inter- 
pret. Figure [3] reveal the terahertz spectra of the elec- 
tromagnon in this geometry. We note at this point that 
similar spectra for this field geometry have been obtained 
previously in Ref. [20|. Similar to the data in Fig. [H 
the spectra in the geometry i?||& reveal an increase of 
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FIG. 4: Comparison of the static and dynamic properties of 
DyMnOa for e\\a and B\\b. Diamonds represent tiie measured 
dielectric contribution of the electromagnon. Circles - dielec- 
tric contribution of the electromagnon as predicted by the 
LST relation. Triangles - static dielectric permittivity with a 
high-frequency value Sao = 25 subtracted. 



the electromagnon intensity in external magnetic fields. 
However, already the comparison of the spectra at 6T and 
at 2T suggests that the increase of the mode intensity is 
not directly correlated with the decrease of the resonance 
frequency. Even without exact analysis of the fits one can 
see that the maxima in £2 for 2T and 6T roughly coin- 
cide in spite of different intensities. Figure H] presents the 
comparison of the static and dynamic properties for the 
geometry B\\b. The electromagnon frequency in this ge- 
ometry softens below 2T and then remains constant (red 
circles). Blue diamonds represent the strength of the 
electromagnon which increases continuously in the whole 
range of the magnetic fields investigated. Contrary to the 
results for the B\\a (Fig. [2]), above the transition to the 
ab cycloid at 2T the mode contribution Ae deviates from 

the LST prediction Ae-^^ (red circles in Fig.©. This 
reflects that the spectral weight of the electromagnon is 
not conserved in external fields i?||6. 

In addition, static dielectric constant as function of 
magnetic field is infiuenced by a second process which 
leads to the peak-like feature as documented in Fig. S) 
Most probably, this feature has to be described using a 
different model. In this case we refer to the recent paper 
[32] by Kagawa et al. in which the dielectric contribution 
of the domain walls in DyMnOa has been investigated. It 
has been shown that the peak in the dielectric constant 
around bc-to-ab phase transition is due to a domain wall 
relaxation with characteristic frequency situated in the 
radiowave range. We conclude that close to i? = 2T 
the main changes in the static permittivity are due to 
the motion of the domain walls. However, except for the 
region close to 2T the overall increase of the static per- 
mittivity corresponds well to the dielectric contribution 
of the electromagnon. As a result, the similarity to the 
soft mode behavior is quite subtle for the B\\b geometry. 
Here a better analogy is probably represented by the mo- 
tion of the domain walls i32l . The relaxation rate of the 



domain wall motion (Fig. 2 of Ref. [321) seem to cor- 
relate well with the peak position of the static dielectric 
permittivity. 



CONCLUSIONS 

The behavior of the electromagnon in DyMnOa in ex- 
ternal magnetic fields has been investigated and com- 
pared to the changes of the static dielectric permittivity. 
A softening of the electromagnon is observed for exter- 
nal magnetic fields parallel to the a-axis, which is ac- 
companied by an increase of the dielectric strength and 
of the static permittivity. Our results demonstrate that 
a soft-mode like behavior with a substantial increase of 
the static dielectric constant can be observed for electro- 
magnons as well. In this case a magnetoelectric mode 
becomes soft and an external magnetic field serves as a 
driving parameter for this behavior. The observed simi- 
larity reveals further close connection between multifer- 
roics and classical ferroelectrics. For magnetic fields par- 
allel to the b-axis no direct correlation between the elec- 
tromagnon frequency and the static dielectric permittiv- 
ity exists. In this case the static properties seem to be 
governed by a motion of the domain walls. 
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